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ABSTRACT 

The low-mass protostellar region L1527 is unusual because it contains observ- 
able abundances of unsaturated carbon-chain molecules including C n H radicals, 
H 2 C n carbenes, cyanopolyynes, and the negative ions C4H~ and CeH~, all of 
which are more associated with cold cores than with protostellar regions. Sakai 
et al. suggested that these molecules are formed in L1527 from the chemical 
precursor methane, which evaporates from the grains during the heat-up of the 
region. With the gas-phase osu.03.2008 network extended to include negative 
ions of the families C~, and C n H~, as well as the newly detected CaN - , we mod- 
eled the chemistry that occurs following methane evaporation at T « 25-30 K. 
We are able to reproduce most of the observed molecular abundances in L1527 at 
a time of ~ 5 x 10 3 yr. At later times, the overall abundance of anions become 
greater than that of electrons, which has an impact on many organic species and 
ions. The anion-to-neutral ratio in our calculation is in good agreement with 
observation for CeH~ but exceeds the observed ratio by more than three orders 
of magnitude for C4H~. In order to explain this difference, further investigation 
is needed on the rate coefficients for electron attachment and other reactions 
regarding anions. 

Subject headings: astrochemistry — ISM: abundances — ISM: molecules — ISM: 
clouds — ISM: protostars: individual(L1527) 



1. Introduction 



L1527, an envelope of a low- mass star-forming region with IRAS04 368+2557 at i ts cen - 
ter, has the physical features of a class 0/1 protostar, as discussed by lAndre et al.l (120001 ) . 
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The infalling collapse of cir cumstellar gas was o bserved by lOhashi et al.l (119971 ). and the 
gas outflow was observed by iTamura et al.l (119961 ) . A number of int erstellar molecules wer e 



discovered in L1527 as part of a molecular survey of dense cores by IjQrgensen et al.l ( 2004). 



Although one might expect the chemistry to be that of a hot corino (jceccarelli et al 
the tem perature of the enve lope is ~ 30 K, well below that of hot corinos (100 K). More 



2003) 



recently, ISakai et al.l (j2008al ) detected unsaturated hydrocarbons and cyanopolyynes in the 
central core region of L1527. These species are unusual in corinos, and are more normally 



assoc ia ted with cold (T 1 K) dark clouds such as TMC-1 (jSuzuki et al.lll992l ; ISmith et al. 



20041 ). ISakai et al.l (j2008al ) explained the production of the carbon-chain species in terms 
of the evaporation of methane during the warm-up to 30 K followed by a rapid gas-phase 
synthetic chemistry, a process that they termed "warm carbon-ch ain chemistry." This chem- 



istry is related to a grain mantle d esorption model for TMC-1 of iMarkwick et al.l (120001 ). In 
a recent paper, lHassel et al.l (120081 ) used the osu gas-grain code to study the warm-up chem- 
istry of L1527 in the vicinity of 30 K and found that the carbon-chain abundances could be 
reasonably reproduced in this "lukewarm corino" by a chemistry in which gas-phase methane 
is the precursor. In addition to these unusual species, the two negative ion species C 4 H~ 
and C 6 H" were also detected in L1527 JSakai et al.ll2007L l2008bl kgundez et al.lboosh . 



An ions i n the interstellar medium were first discussed by iDalgarno fc McCrayl (119731 ) . 



Herbstl (Il98ll ) proposed that neutral radicals with large enough electron affinities can undergo 
efficient radi ative attachment to form anions in cold dense clouds if they possess more than 
« 4 atoms. iLepp &: Dalgarnol (119881 ) suggested that negatively c harged polycyclic aromatic 



hydrocarbons (PAHs) could also form by radiative attachment. IWakelam fc Herbstl (120081 ) 
modeled the gas-phase chemistry in dense cold cores with neutral and negatively charged 
PAHs, and found that the inclusion of PAH's allows good agreement with observations for 
small species with the use of so-called "high-metal" elemental abundances. Confirmation of 
anions in interstellar and circumstellar sources awaited the measurement of their rotational 
spect ra in the laboratory by Thaddeus and co-workers ( McCarthy et al.l 20061; Gupta et al. 
2007h. In addition to the detections in L1527, C 6 H" and C«H~ JMcCarthy et al.l bood: 
Brunken et aD l20Q7h were detected in TMC-1, while dH" JCernicharo et alT bo07h. C 6 H" 



(IMcCarthy et al.l 120061 ; iKasai et al.l 120071 ) . and C 8 H (IRemijan et al.l 120071 ) were observed 



and C„,H 



n=4- 



in IRC+10216. In response to these detections, the anions C n , n=5-10, 
10, w ere added to the RAT E06 (udfa.net) and IRC+10216 chemical net works dMillar et al. 



2000) by lMillar et al.l (120071 ) with estimated or calculated rate coefficients ( iTerzieva fc Herbst 



20001 ) for radiative attachment and photodetachment processes as well as measured rate 
coefficients for associative detachment reactions. The augmented models were then used to 
study anionic abund ances in TMC-1, IR C+10216, and photon-dominated reg ions suc h as th e 
Horsehead Nebula. iMillar et al.l (120071 ) found that, as predicted earlier by iHerbstl (Il98ll ). 



- 3- 



hydrocarbon radical anions can reach abundances exceeding 1% of their neutral precursors, 
depending on their size. Another type of a nion, C3N - , the radiative attachment of which 
was first studied by IPetrie fc Herbstl (119971 ). was also detected recently in IRC+10216 by 



Thaddeus et alj (2008). Its radiative attachment rate coefficie nt has been calculated anew 



by iHerbst fc Osamural (120081 ) , with a result similar to that of IPetrie fc Herbstl (119971 ) used 
here. 

In this paper, we re-study the gas-phase chemistry of L1527 with particular interest in 
the formation and depletion of anions and how their chemistry affects the overall chemistry of 



the lukewarm corino. The osu.03.2008 gas-phase network (see http:/ /www. physics. ohio-state.edu/ ~eric/ ) 
has been augmented to include formation and depletion reactions for the carbon cluster and 
hydrocarbon radical anions as well as for C3N - . The remainder of this paper is organized as 
follows. In Section [21 we discuss the gas-phase code used and the choice of initial abundances. 
In Section [3j the results are discussed and compared with observation for both anions and 
other species and with the earlier gas-grain results for non-anionic species. In Section HI we 
consider how our results respond to variations in initial abundances and to the exclusion of 
anionic chemistry. Finally, we end with a summary. 



2. Gas-Phase Code and Initial Conditions 



The recent osu gas-phase network (http:/ /www. physics. ohio-state.edu/~eric/ ) osu.03.2008 
has been augmented by the addition of several families of negative ions and classes of reac- 
tions to form and destroy them. The new negative ions included in the model consist of the 
carbon chain families C~ (n = 5-10) and C n H~ (n = 4-10) and the newly detected anion 
C3N - . The reaction classes include radiative attachment (A + e — > A~ + hi/), photode- 
tachment (A - + hi/ — > A + e), dissociative attachment (AB + e — > A~ + B), associative 
detachment (A - + B — ► AB + e), anion-neutral reactions (A~ + B — > C~ + D) and 
mutual neutralization (A~ + B + — > A + B). Most of these class es of reactions are con- 
tained in the modified RATE06 network used in iMillar et al.l (120071 ). We have updated the 
rate coefficients for the r adiative attachment of C4H and C5H according to the theory of 
Herbst fc Osamural (120081 ) . A large number of mutual neutralization reactions involving all 
anions an d 57 major positive ions were added with rat e coefficients k = 7.5 x 10 _8 (T/300) -0 5 
cm 3 s -1 ( ISmith et al.lll978l ; IWakelam fc Herbstl 120081 ); without these reactions, the correct 
balance of anions, electrons and positive ions cannot be maintained for L1527 at later times, 
when the overall anionic abundance exceeds the electron abundance. Newly added reactions 
for the formation and depletion of C3N - and their rate coefficients are listed in Table [TJ 
Overall, the augmented network for L1527 contains 5941 reactions involving 469 species, an 
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extension of osu.03.2008, which contains 4479 reactions involving 468 species. 



We follow the gas-phase chemistry after its warm-up from a cold cloud with icy mantles 
to a temperature in the range 25 — 30 K, by which most of the methane ice in grain mantles 
has evaporated along with other volatile species such as CO and N2, which evaporate at 
slightly lower temperatures than methane. We u se a temperature of 30 K, a total densit y 
nn = 10 6 cm -3 and visual extinction A v =10 mag (ISakai et al.ll2008al ; iJprgensen et al.ll2004l ). 
To convert observed co lumn densities into fra ctional abundances, we use a column density 
Nh = 6 x 10 22 cm -2 (iJprgensen et al.l 120041 ) . The cosmic ray ionization rate ( is taken 
to be 1.3 x 10 



-17 



The initia 



abu ndances, shown in Table [2J are based on the gas- 



grain code results by lHassel et al.l (120081 ) . in which the gas and grain populations during the 



warm up are calculated. Since most of the water ice resides on grain mantles until much 
higher temperatures, the elemental carbon-to-oxygen abundance ratio in the gas is carbon 
rich, as can be seen by adding up the carbon-containing and oxygen- containing species. The 
important precursor methane has a fractional a bundance of 3 x 10 ~ 6 , in reasonab le agreement 



with values obtained for cold ices using ISO ( jGibb et al.l 120041 ) and Spitzer ( jOberg et al. 



20081 ). The low abundances for metals reflect the fact that most are still embedded in grain 
mantles. The chemistry starts quickly, as methane is pa rtially converted in to acetylene and 
more complex hydrocarbons via ion- molecule reactions (IHassel et al.l 120081 ). 



3. Results 

3.1. Hydrocarbon Radicals and Anions 

Hydrocarbon radicals of the general linear structure C n H are the main precursors to 
the C n H~ anions through radiative attachment. The calculated fractional abundances (with 
respect to n H ) of C 4 H, C 6 H, C 8 H, and Ci H along with their respective anions as functions 
of time are shown in Figure [TJ For the first 3 x 10 4 yr, the neutral radicals increase in 
abundance with abundances inversely proportional to size. At 10 3 yr, the main production 
channels of C n H are dissociative recombination from precursor ions C n H+: 

C n H++e" — C n H+..., (1) 
and neutral-neutral reactions from smaller hydrocarbons: 

C + C n -iH 2 — > C n H + H. (2) 
The main formation reactions after 10 4 yr become the associative detachment processes 

C- + H — > C n H + e" (3) 
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because the concentrations of the bare cluster anions become significant. When peak values 
are achieved, at a time of ~ 3 — 5 x 10 5 yr, the size vs fractional abundance order has 
broken down, and the most abundant radical becomes C§H. The radical CgH is next, while 
the still lower abundances of C4H and C10H are similar. As steady state is achieved at still 
later times, the abundances of the three smallest radicals are nearly the same, while C10H is 
somewhat lower. The unusual dependence on size is due partially to carbon-rich elemental 
abundances and partially to the inclusion of anions, which increase the synthetic power of 
the chemistry and cause an "edge" effect, as discussed below in Section 1472"! 

The calculated fractional abundances of the C n H anions, depicted in Figure [U, show an 
unusual dependence on size as well. At times after 4 x 10 4 yr, the most abundant anion is 
CgH~, which slightly exceeds the smaller CqH~. There is a gap in abundance between these 
two, and the pair of ions CioH - and C4H - , which trail by a factor of 3-5. The relatively 
low conc entration of C/iH~ is cau sed by the relative slowness of the radiative attachment 
reaction (IHerbst fc Osamurall2008l ) 



C 4 H + e" 



C 4 H" + hi/, 



(4) 



despite the fact that the calculated value of the rate coefficient is known to be too large 
( IHerbst fc Osamurall2008l ). The observed frac tional abundances ofC^H - and CeH - in L1527 
are 1.8 x 1(T 13 and 9.7 x 1(T 13 respectively JSakai et alil2007L l2008ah . It can be seen that 
while the calculated value for the larger anion is close to the observed value at relatively 
early times (~ 10 3-4 yr), the calculated abundance for the smaller anion is far too large at 
all but the shortest times. 



3.2. Cyanopolyynes and C3N 

Cyanopolyynes are produced via a variety of ion-molecule and neutral-neutral channels. 
The existence of the C n H~ anions boosts the efficiency of one production mechanism: 

CN + C 2n H 2 — >HC 2n+1 N + H (5) 

by enhancing the abundance of the the precursor C211H2 through the associative detachment 
reactions 

C 2n H~~ + H — > C 2 nH 2 + e~. (6) 

Although cyanoacetylene, HCCCN, can be formed by neutral-neutral channels, its higher 
energy is omers, such as HNCC C, are thought to be formed only through dissociative recom- 
bination (jOsamura et al.lll999l ); viz., 

C 3 H 2 N+ + e" — ► HNCCC + H. (7) 
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This particular isomer is important, because as can be seen in Table [U, it leads to the 
anion CaN - via an exothermic dissociative attachment reaction with hydrogen atoms. This 
synthesis can be competitive with the more simple radiative attachment depending on the 
source. Indeed, in L1527, it is the dominant pathway, as indicated in Figure [2], which shows 
the calculated fractional abundances of HC3N, HNCCC, C3N and CsN~ as functions of time. 
It can be seen that the abundance of the radical C3N exceeds that of the isomer HNCCC for 
all times past 2xl0 3 yr. However, the abundance of the radical does not exceed that of the 
isomer by more than a factor of 3, while the radiative attachment rate coefficient for C3N - 
+ e~ is lower than the dissociative attachment coefficient for HNCCC + e~ by a factor of 
100, as shown in Table [TJ The predicted fractional abundance of the anion lies somewhere 
in the range 10~ n — 10~ 12 for most times. Although this range of values seems very low, it 
should be noted that the anion C4H~ was detected with an even lower fractional abundance. 



3.3. Anion-to-Neutral Ratio 



The calculated and observed anion-to-neutral abundance ratios vs time are shown in 
Figure[3]for C4H and CqB.. Despite the fact that both anion and neutral are significantly over- 
produced at times > 10 4 yr, the calculated C 6 H~/C 6 H ratio always lies close to the observed 
value of ~ 10%, exceeding it by at most a factor of five at early times, and dipping below 
by less than a factor of two before steady state is reached. This level of agreement suggests 
that the theoretical rate coefficient for the radiative attachment of CgH ( iHerbst fc Osamura 
20081 ) is reasonably accurate although the alternative possible formatio n route via disso- 
ciative attachment of th e carbene H 2 Cg is not included in our network (ISakai et al.l 120071 ; 
Herbst fc Osamural 120081 ) . The ratio C 4 H~/C 4 H, as for other sources, is calculated to be 
much too large, due principally to the failure of the phase space treatment for the radiative 
attachment of C4H ( Herbst &: Osamural 2008 ). 



In Figure HJ we plot the fractional abundances of the electron, total positive ions, and 
total anions. After 4 x 10 4 yr, the total anion abundance exceeds the electron abundance, 
and anions become the dominant form of negatively charged particles. Since the mutual 
neutralizations of positive ions and negative ions are less rapid than dissociative recombi- 
nation reactions between positive molecular ions and electrons, negative ions are likely to 
stay more abundant than electrons. Indeed, at steady state, the anionic abundance exceeds 
the electronic abundance by a factor of more than two. Since anions can be more abundant 
than electrons, it was imperative to add a large number of mutual neutralization reactions 
to make sure that the correct proportion of anion, electron, and positive ions was obtained. 
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3.4. Overall Comparison with Observation 

To determine the best fit to the observed 23 abun dances for the anions and other species 



in L1527, we used the mean confidence level method of iGarrod et al.l (120071 ). In this method 



a confidence level Ki is defined for each species i by the equation 

log(Xi) - log(X obSji ) 



^ = erfc 



V2< 



a 



where X, is the computed fractional abundance, X b s ,i the observed fractional abundance, 
and (7 = 1. The agreement is perfect for factor of three departure from observation 

if Ki = 0.63, and a one order-of-magnitude departure from observation if Ki = 0.32. We then 
take k, the average value of Ki, as the criterion for agreement. In our network of reactions, 
we do not make a distinction between the acetylenic forms HC n H and the carbenes H2CV 
Since only the latter are d etected, we make the assumption that they represent only 1% of 



the total C n H 2 abundance (IPark et al.ll2006l ). The optimal value of k is 0.652 corresponding 
to an average agreement within a factor of three; this agreement occurs at a time of 4.8 x 10 3 
yr. The species and their observed and optimal calculated fractional abundance are shown in 
Table [31 where it can be seen that 18 of the 23 molecules have computed abundances within 
an order of magnitude of the observed values, while 15 have computed abundances within a 
factor of three of the observed values. The most extreme outlier is the anion C4H~, which is 
predicted to have an abundance higher than observed by a factor of 61. The other molecules 
that disagree with observation by more than an order of magnitude are C4H, N 2 H + , and 
CN. Another anion, CeH - is overproduced by an order of magnitude, C5H and HNC are 
overproduced by more than a factor of a few, while HCO^ is underproduced by factor of 
five. The mean confidence limit is not a strong function of time despite the fact that the 
abundances of carbon-chain molecules tend to keep increasing past the optimal time; at 
steady-state, which occurs at times slightly in excess of 10 6 yr, the value of k is still 0.474, 
corresponding to an average discrepancy factor of slightly more than five. At steady-state, 13 
species are in agreement within an order of magnitude, and 9 of them are in agreement within 
a factor of 3. The main difference between the optimal time and steady-state results concerns 
the abundances of carbon chain species, which become much larger than their observational 
values past the optimal time. In fact, at 5xl0 5 yr, when these molecules achieve their peak 
abundances, the confidence level is 0.411; a little worse than at steady-state. 

If we exclude the anions from our network, and fit for the remaining 21 molecular 
abundances in L1527, we find an optimal value for k of 0.628 at 5.9 x 10 3 yr, corresponding 
to an average factor of three discrepancy, and a steady-state value of 0.487, corresponding to 
an average factor of five discrepancy. The optimal results without anions are a lso tabulated 
in Table [3] as are the optimal results of the gas-grain warm-up model to 30 K (jrlassel et al. 
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20081 ). This latter model, which does not contain anions, achieves an optimal agreement of a 
factor of three discrepancy when the temperature has risen to 25 K from its initial value of 
10 K. The time spent up to 25 K when the gas-phase carbon chemistry can occur is roughly 
1 — 2 x 10 4 yr. Comparison of the gas-phase models with the gas-grain model shows that 
the gas-phase model without anions achieves somewhat closer agreement with the gas-grain 
model, because, as discussed later, the anions do affect the chemistry. It should also be 
noted that the observed molecular abundances towards LI 52 7 excluding anions can also be 
fit to within an order of magnitude by a cold model, especially the non-c arbon-chain species , 



which may not be spatially correlated with the carbon-chain molecules (IHassel et al. 



20081 ). 



Predicted abundances for a wide class of molecules using the optimal times of the gas- 
phase fit with anions as well as the gas-grain warm-up model to 30 K are listed in Table HI 
Most optimum gas-phase abundances tend to be only slightly smaller than their gas-grain 
counterparts. Other than the negative ions, some significant disagreements (> 1 order-of- 
magnitude) concern the hydrocarbon radicals, which are depleted by electron attachment in 
the gas-phase model, the carbenes, and species such as methanol, which are not produced 
efficiently in the gas phase. Despite the fact that the gas-grain prediction for methanol is 
much higher than the gas-phase result, methanol is still predict ed to be of very low abundance 



in the gas until the temperature warms up to hot core levels (IHassel et al.ll2008l ). 



4. Response to Variations 



4.1. Varying Initial Conditions - Oxygen and Methane 



Unlike the gas-grain approach of IHassel et al.l (120081 ) to L1527, where the evaporation 
of the volatile mantle species is determined by the chemical and physical processes, in the 
gas-phase approach we can vary the initial abundances from those listed in Table [2J which 
derive from evaporation from grain mantles in the warm-up gas-grain model. It is useful to 
be able to vary the initial abundances because the warm-up scenario is approximate at best. 
Here we consider variations in the initial abundances of atomic oxygen and methane (CH4). 
The C/O elemental abundance ratio is 1.06 with our standard abundances; an increase in O 
will move the ratio closer to unity. The amount of methane is important to the warm carbon 
chemistry, since it is the pre cursor for all of it, and since its observed prior abundance in ices 
can vary (jOberg et al.ll2008l ). 



Table [5] shows the effect of changing the initial abundance of atomic oxygen on the 
time of best agreement, the average confidence level k, and the number of species within 
order-of-magnitude agreement with observation. Although the O abundance is not large, O 
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atoms tend to destroy organic ions and radicals. Varying the fractional O abundance from 
5 x 10~ 8 to 1 x 10 -6 shows that the optimal value for k peaks for 1 — 2 x 10~ 7 , although the 
criterion of number of molecules fitted to within an order-of-magnitude peaks at 5 x 10~ 7 . 
The sensitivity to changes in O abundance, however, is minimal. On the other hand, the time 
of best agreement increases with increasing initial abundance of O. This latter dependence 
arises from the fact that O slows down the synthetic process. 

Regarding methane, the results for time of best agreement, k, and order-of-magnitude 
agreement are shown in Table There is little dependence on the optimal fit defined by 
either the k criterion or the number-of-molecules criterion as the initial methane abundance 
is raised from 3 x 10 -7 to 3 x 10 -5 . Here, though, the time of optimal agreement decreases 
with increasing methane abundance. Presumably the higher abundance of methane leads to 
the more rapid production of carbon-chain species. 



4.2. The Effect of Anions 



It is interesting to compare the chemistry of our gas-phase model for L1527 with and 
without anions in more detail. Given that the overall abundance of anions in the model 
exceeds that of electrons at times later than 4 x 10 4 yr, one would expect the effect of anions 
to be more pronounced at such times. By and large, this inference is true although for 
some classes of molecules there is little effect at all at earlier times, while for other species 
the effect is already a large one. In Figure [5j we plot the fractional abundances of four 
classes of molecules vs time - the C 2n H radicals (Upper left panel), the hydrocarbons C n H 2 , 
consisting mainly of polyynes (HC 2n H) with small abundances of carbenes H 2 C 2n (Upper 
right panel), the bare carbon clusters (Lower left panel), and the cyanopolyynes HC 2n +iN ( 
Lower right panel) - with and without anionic chemistry. Let us start with the hydrocarbons. 
Here it can be seen that the effect caused by the anions increases with time, as expected, 
so that by steady-state the abundances computed without anions are mainly significantly 
lower than those computed with anions. The reason is that the C n H anions increase the 
abundances of the hydrocarbons via reaction (J6j) as well as the analogous processes for odd 
numbers of carbon atoms. We also note that at later times, including steady state, the largest 
hydrocarbon depicted, CsH 2 , has the largest abundance in the model with anions, and the 
smallest abundance in the model without anions. The chemistry leading to this unusual 
result is complex and will be the subject of a future paper. For now, we note only that the 
existence of anions has a catalytic effect on the production of more complex hydrocarbons 
from simpler ones. For example, starting with the hydrocarbon CeH 2 , reaction with carbon 
atoms leads to the species C7H (see reaction (jSJ)), which is efficiently converted to its anion. 
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The anion, in turn, undergoes an associative detachment reaction with atomic hydrogen to 
form the more complex hydrocarbon C7H2 (see reaction (j3J)). The process continues through 
the largest hydrocarbon in the model. The enhanced abundance of this largest hydrocarbon 
suggests that still larger species must be included in the model to avoid an "edge" effect 
for potentially observable species, in which the concentrations of the larger species are over- 
estimated due to a paucity of destruction reactions. 

If we look at the cyanopolyynes, the result, though more complex, is similar to that 
of the hydrocarbons. The enhanced synthetic power with anions is doubtless caused by 
the enhanced production of these hydrocarbons, which can react with the radical CN to 
produce cyanopolyynes, as shown in reaction (jSJ). The late-time abundance of the largest 
cyanopolyyne shown, HC13N, is comparable with those of smaller cyanopolyynes. On the 
other hand, the plot of the C^nH radicals shows that the existence of anions decreases their 
abundances starting at rather early times. This effect is surprising because bare carbon an- 
ions are eventually involved in their formation (see reaction ([3])), but presumably destruction 
via the formation of anions is more important. For the bare carbon clusters, destruction by 
the formation of anions causes a similar effect. 



5. Summary 



LI 527 is an unusual source because the unsaturated carbon-chain molecules detected 



there a re more normally associated with cold dense cores at lower temperatures. ISakai et al. 



( 12008al ) proposed that the carbon-chain chemistry starts soon after rising temperatures allow 
the volatile species methane to evaporate int o the gas. Ou r grou p has studied this chemistry 
via two types of models. In a previous model, lHassel et al.1 (120081 ) used the osu gas-grain code 
to study the chemistry in a cold dense core followed by warm-up to temperatures associated 
with L1527 (T ^ 30 K) and found a rich chemistry leading to carbon-chain species, as 
hypothesized by ISakai et al.1 (j2008af ). Comparison with observation showed that although 
the non-carbon-chain species in L1527 could be fit well either as remnants of a previous cold 
phase or current denizens of the warm-up phase, the carbon-chain species are fit significantly 
better with the warm model. Although the gas-grain approach was quite useful, this method 
does not contain any anionic chemistry, and the chemistry is fixed by what happens in the 
earlier cold epoch and the rising temperatures. We have followed this prior study with the 
isothermal (T = 30 K) gas-phase study reported here because it allows us to include the 
gas-phase anionic chemistry without worrying about what happens when anions collide with 
grains and because it allows us to vary some of the initial abundances, which are mainly 
based on the earlier model. The abundance of methane, for example, has been studied in 
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the ice phase of cold sources (lOberg et al.ll2008l ) prior to evaporation and found to depend 
somewhat on the source investigated. 

Including the chemistry of carbon-chain anions in our gas-phase network has led to a 
large increase in the number of reactions, mainly because of the many cation-anion neu- 
tralization reactions needed. With the chosen initial abundances, the best agreement with 
all 23 observed molecules in L1527 occurs at a time of 4.8 x 10 3 yr, when the abundances 
of most carbon-chain species are increasing. This agreement corresponds to an average 
factor of three discrepancy between observational and calculated values. Variation in the 
initial abundances of either atomic oxgyen or methane changes the time of best agreement 
mainly because these species regulate both the formation and depletion rates of carbon-chain 
species. For the observed anions, C4H~ and CgH~, we fit both the fractional abundance and 
the anion-to-neutral ratio reasonably well for the larger anion, but dramatically overpro- 
duce these quantities for the smaller anion. Finally, it is to be noted that the steady-state 
abundances in the gas-phase model are reached by 10 6 yr, so that the chemical processes of 
interest in this paper are completed well before the gas-grain models show declines in the 
gas-phase abundances (iHassel et al.ll2008l ). 



A comparison of the chemistry occurring with and without negative ions shows that 
anions have a major effect on the chemistry of other species, especially at later times than 
that for the optimal fit, when the overall anionic abundance exceeds that of electrons. This 
effect can either lower abundances of a given class of molecules or raise them, depending upon 
whether the formation of anions tends to deplete other species or enhance their production. 
There is also a strange effect in which the dependence of the abundance of a molecule in 
a family such as the cyanopolyynes on size is changed when anions are included in the 
chemistry. Typically, in the absence of anions, an increase in size results in a drop in peak 
abundance. The inclusion of anions appears to make this dependence significantly smaller 
and even to vitiate it completely in some instances. As will be shown in subsequent work, 
this effect is also pronounced in the cold core TMC-1 under oxygen-rich conditions. With 
relatively large abundances for the largest carbon-containing species in our network, it may 
well be necessary to add still larger molecules to the model, either by including all of the 
formation and depl etion reactions that govern their chemistry, a s in the IRC+ 10216 network 
(IMillar et al.ll2000l ). or by some more approximate approach (IHerbstl Il99ll ). In addition, 
despite the size of the current network, the role of negative ion- neutral reactions remains to 
be fully explored JCordiner et alil2008l : ICordiner k Millar! l2008f ). 
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Table 1: Formation and Depletion Reactions for C3N 
Reaction Rate Coefficient (cm 3 s _1 ) Reference 



C 3 N + e -> C 3 N" + hi/ 2.0 x 10- 10 (T/300)~ a5 1 

HNCCC + e -> C 3 N - + H 2.0 x 10- 8 (T/300)-°- 5 2 

C 3 N~ + H -f HC 3 N + e 1.0 x 10" 9 3 

C 3 N~ + C -> C 4 N + e 1.0 x 10- 9 3 

C 3 N~ + O -> C 2 N + CO + e 1.0 x 10" 9 3 

C 3 N" + hp -> C 3 N + e 2.0 x 10~ 9 exp (-2.0A V ) 4 

C 3 N~ + C+ -> C 3 N + C 7.5 x 10- 8 (T/300)-°- 5 5 



References. — (1) IPetrie fc Herbsti (119971) 



(2) estimate based on 



Adams et al 



(198' 



(3) based on Eichelberger et al (|2007f) 

(4) based on lMillar et al.l hoO^I 



(5) sample positive ion; see lSmith et alJ (|1978l ) for typical rates 
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Table 2. Initial Abundances 



Species 


Fractional Abundance 


Species 


Fractional Abundance 


C 


1.0(-08) 


He 


6.0(-02) 





1.0(-07) 


Fe 


l.O(-ll) 


Na 


l.O(-ll) 


S 


1.0(-08) 


Si 


3.0(-13) 


CO 


5.0(-05) 


cs 


3.0(-09) 


CH 4 


3.0(-06) 


H 2 


5.0(-01) 


CN 


4.0(-09) 


HCN 


3.0(-09) 


HNC 


3.0(-09) 


H 2 


1.0(-08) 


H 2 S 


1.0(-09) 


N 2 


1.0(-05) 


NH 3 


1.0(-08) 


o 2 


1.0(-08) 


C+ 


1.0(-09) 


He+ 


l.O(-ll) 


H 3 + 


6.0(-ll) 


S+ 


l.O(-ll) 


Fe+ 


l.O(-ll) 


Na+ 


l.O(-ll) 


HCO+ 


1.9(-09) 


e~ 


3.0(-09) 


(C)elem 


5.3(-05) 


(O)elem 


5.0(-05) 







Note. — a(-b) represents ax 10 b . 
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Table 3. Observed and Optimal Calculated Abundances in L1527 



Species 


Observed Abundance 


Reference 


Calculated Abundance 


Without anions 


Gas-Grain Model 


CO 


3.9(-05) 


(1) 


5.0(-5) 


5.0(-5) 


5.1(-05) 


CN 


8.0(-ll) 


(1) 


1.4(-9) 


l.l(-9) 


9.6(-10) 


CS 


3.3(-10) 


(1) 


5.5(-10) 


5.9(-10) 


5.1(-09) 


C 2 S 


8.5(-ll) 


(2) 


1.6(-10) 


1.7(-10) 


1.3(-10) 


HCN 


1.2(-9) 


(1) 


1.3(-9) 


1.3(-9) 


7.7(-10) 


HNC 


3.2(-10) 


(1) 


1.5(-9) 


1.4(-9) 


6.3(-10) 


SO 


1.4(-10) 


(1) 


2.7(-10) 


3.3(-10) 


1.3(-10) 


C 4 H 


3.3(-9) 


(2) 


1.3(-10) 


5.7(-10) 


1.7(-09) 


C 5 H 


1.6(-11) 


(2) 


9.5(-ll) 


1.7(-10) 


4.6(-10) 


C 6 H 


l.O(-ll) 


(2) 


1.9(-11) 


4.8(-ll) 


2.3(-10) 


C 4 H- 


1.8(-13) 


(3) 


l.l(-ll) 






C 6 H" 


9.7(-13) 


(4) 


9.8(-12) 






HC 3 N 


8.9(-10) 


(1) 


5.5(-10) 


5.8(-10) 


4.0(-10) 


HC 5 N 


9.7(-ll) 


(2) 


l.l(-lO) 


1.4(-10) 


1.6(-10) 


HC 7 N 


2.7(-ll) 


(2) 


1.4(-11) 


1.8(-11) 


2.9(-ll) 


HCgN 


2.5(-12) 


(2) 


2.1(-12) 


2.3(-12) 


5.9(-12) 


H2C3 


l.O(-ll) 


(2) 


7.9(-12) 


8.7(-12) 


l.l(-ll) 


H2C4 


2.2(-ll) 


(2,5) 


1.7(-11) 


2.2(-ll) 


1.2(-10) 


H2C6 


2.5(-12) 


(2,5) 


3.0(-12) 


3.7(-12) 


2.5(-ll) 


C3H4 


7.8(-10) 


(2) 


1.4(-9) 


1.7(-9) 


2.2(-09) 


HCO+ 


6.0(-10) 


(1) 


1.3(-9) 


1.2(-9) 


1.7(-09) 


N 2 H+ 


2.5(-10) 


(1) 


l.O(-ll) 


l.O(-ll) 


l.l(-ll) 


HCO+ 


1.0(-12) 


(6) 


1.9(-13) 


1.8(-13) 


4.1(-13) 



^J^rgensen et alj l|2004h 
feakai et al] ll2008ah 
^Agundez et all ll2008ri 
iSakai et all j2008bh 

5 Calculated carbene abundance was adjusted to 1% of total C n H2 (see Section 3.4). 
ISakai et all J2008ch 



Note. — a(-b) represents ax 10 . 
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Table 4. Predicted Optimal Fractional Abundances of Unobserved Gaseous Species 



Species Gat 


5-Phase Model 


Gas-Grain Model 


p w 
U 3 ri 


o.oy-iv ) 


±.u^-y ) 


P^H 


A A( 1 91 


1 3f inl 

1 . 1 1U J 




AM 1 91 


1 1 ( 1 nl 


G„H 


8 "\( 1 31 

O.OI id J 


7 9f 111 


CioH 


4 *\( 1 31 
Hi.oy lo j 






3 f\( 111 

O . U I X J- J 






7 9f 1 91 






1 Q( 1 9^1 










P 


O (1/ 1 Ql 




PTT 


/.y ^-0 j 


o.y ^-0 j 


C2H2 


1 4f Si 


3 7f 81 


H2C5 


fi n/ 1 91 


9 ^( 111 

Z.O\ - L L 1 


HoG^ 
£12^7 


1 4f-19l 


7 1( 1 91 


H2C8 


^ 1 31 


^ fif 1 91 

O.UI li 1 


H2C9 


9 fif 1 31 


9 fir 1 91 

i.UI li J 




9 nr 1 ni 


1 1 ( Ql 




1 A( 111 
1.41-11 1 


1 7f inl 




1 1 f 1 31 

J- . -L I 1U 1 






O . Z\ LO 1 




P 1\T 


0/ 1 n\ 
z.zy- ru J 


/i if ~\r\\ 
4.0^-iu ) 


PoN 


fi nr 1 1 1 


1 9( 1 nl 

l.i^-lU J 


C 4 N 


4.3(-ll) 


1.7(-10) 


C 5 N 


8.6(-ll) 


1.4(-10) 


C 7 N 


l.O(-ll) 


2.6(-ll) 


C 9 N 


1.3(-12) 


5.1(-12) 


CH3C3N 


1.5(-13) 


2.9(-13) 


CH 3 C 5 N 


7.5(-15) 


2.3(-14) 


CH 3 C 7 N 


8.2(-16) 


3.8(-15) 


C 3 N" 


5.5(-12) 




N 2 


1.0(-5) 


5.6(-6) 
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Table 4 — Continued 



Species 


Gas-Phase Model 


Gas-Grain Model 


o 2 


7.4(-9) 


1.1 (-8) 


H 2 


7.4(-10) 


8.8(-10) 


NH 3 


7.0(-9) 


3.7(-9) 


co 2 


2.2(-9) 


4.1(-9) 


He+ 


1.9(-12) 


8.8(-12) 


C+ 


2.9(-ll) 


l.l(-lO) 


HCO+ 


1.3(-9) 


1.7(-9) 


H3 


2.7(-ll) 


5.6(-ll) 


h 0+ 


8 8/ 19^ 
o.o^-lZ J 


1.0^-11 ) 


c-C 3 H 2 


4.2(-10) 


5.7(-10) 


CeH 6 


6.5(-ll) 


3.2(-10) 


CH 2 CN 


4.4(-ll) 


7.4(-ll) 


H 2 CO 


5.4(-9) 


4.3(-9) 


HCOOH 


3.1(-12) 


3.6(-12) 


CH 2 CO 


1.2(-9) 


1.0(-9) 


CH 3 CN 


2.8(-12) 


2.0(-12) 


CH3OH 


3.4(-16) 


5.0(-ll) 


HCOOCH3 


1.8(-19) 


3.8(-15) 


CH3OCH3 


2.3(-17) 


3.8(-16) 


H 2 CS 


1.2(-9) 


5.2(-10) 


H 2 S 


7.7(-10) 


3.3(-9) 



Note. — a(-b) represents ax 10 b . 
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Table 5. Dependence of Agreement with Observation on the Atomic Oxygen Abundance 



Initial fractional 


K 


Time of 


Number of Species 


Oxygen Abundance 




Best Agreement (yr) 


in Agreement 


5.0(-08) 


0.649 


4.3(03) 


18 


1.0(-07) 


0.652 


4.8(03) 


18 


2.0(-07) 


0.652 


5.9(03) 


19 


5.0(-07) 


0.636 


7.2(03) 


20 


1.0(-06) 


0.621 


8.0(03) 


18 



Note. — a(-b) represents ax 10 b . 



Table 6. Dependence of Agreement with Observation on the Methane Abundance 



Initial fractional 


K 


Time of 


Number of Species 


Methane Abundance 




Best Agreement (yr) 


in Agreement 


3.0(-07) 


0.601 


8.8(03) 


18 


1.0(-06) 


0.642 


6.5(03) 


18 


3.0(-06) 


0.652 


4.8(03) 


18 


1.0(-05) 


0.652 


4.3(03) 


19 


3.0(-05) 


0.647 


4.3(03) 


19 



Note. — a(-b) represents ax 10 
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Fig. 1. — The fractional abundances with respect to hh for C211H and C2 n H (n = 2-5) are 
plotted vs time for our model of LI 527. 



-22 - 




Fig. 2. — The calculated fractional abundances with respect to n H for C 3 N , HCCCN, 
HNCCC, and C 3 N are plotted against time for L1527. 
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Fig. 3. — The calculated anion-to-neutral abundance ratio for C 4 H and C 6 H as a function 
of time. The observed values are plotted as horizontal lines. 
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Fig. 4. — The calculated fractional abundances of electrons, the total positive ions, and the 
total negative ions are plotted against time. Note that the total negative ion abundance 
exceeds the electron abundance from a time of 3 x 10 4 yr. 
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Fig. 5. — Comparison of the time-dependent abundances of selected families of molecules - 
hydrocarbon radicals (C2 n H; Panel 1), hydrocarbons (C211H2; Panel 2), carbon clusters (C2 n ; 
Panel 3), and cyanopolyynes (HC 2n +iN; Panel 4) - with and without negative ions. Solid 
lines show the results with negative ions, while dotted lines show the results without negative 
ions. Note that CiqH is not included in the network without anions. 



